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1. INTRODUCTION 


This report describes the computerized procedures which have been 
developed for the protection coordination of the Kennedy Space Center 
Electric Distribution Network. The developed procedure enables the user 
to: 

(1) Visualize the coordination and protection quality of any 
existing system of devices and settings by computer plotting 
the tripping characteristics of the involved devices on a 
common basis, 

(2) Determine the optimum settings of a given set of protective 
devices and configuration in the sense of the best expected 
coordinated operation of these devices. 

In the first case the computer procedure aqcepts data on the network 
configuration, protective devices and settings and provides a computerized 
plot of the overlayed current-time characteristics on a common scale. 

Thus the tedious hand plotting procedure is eliminated* For the KSC 
system, all basic protective elements have been computer modeled in the 
form of subroutines so that input device data is minimal - being only 
settings for adjustable elements. 

The second procedure provides an optimization feature which repults 
in a best coordination available from candidate systems. The user speci- 
fies the candidate devices and available settings of these devices. The 
computer program then compares all combinations of specified equipment 
with regard to coordination and selects the best set together with 
adjustment settings. If only one set of devices is specified, the pro- 
gram selects optimum adjustment settings for this set. The program 
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utilizas the Short Circuit Program as a subroutine to determine network 
fault currents to input to the device characteristic subroutines. 

The developed procedure consists of a set of computer subroutines 
referred to as the coordination pack. The coordination pack is divided 
into five (5) sets of computer routines 
(a) DEVICE cmMCTERISTIC SUBROUTINES 


These are subroutines which simulate the time versus current 


char act erist ic s 

of the different protective devices existing 

in the system. 

The following subroutines are included in this 

set. 


SELIA.C 

- relays type IA.C 

EEIAGO 

- relays type CO 

BRKEAK 

•* breakers t 3 rpe AK 

BRKSIA 

- breakers type lA 

BRKROD 

- breakers type OD 

BRKRDB 

- breakers type DB 

MGCBGE 

- molded case breakers by 6,E. 

FOSEEJ 

- fuses type EJO, EJ 

EUSGLE 

- fuses t 3 ^e CLE 

FUSEGL 

- fuses type GL 

FUSESM 

~ fuses type SM 

FDNCTN 

- this subroutine simulates a 
piecewise first degree function, 
and is used by the above sub- 
routines 

The above 

subroutines return the time to trip required 


by the corresponding device when the current,, is specified. 
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The user who becomes acquainted with these subroutines and I 

with their format may obtain a computer representation of ; 

i 

the particular device by calling the appropriate subroutine. | 

(b) GOORDIHATION INDEX COMPUTATION SUBROUTINES I 

This set of subroutines includes those which calculate | 

i 

probabilistic indices of coordination; namely, it includes j 

I 

the following subroutines; ^ 

\ 

FUSBEK - calculates coordination indices j 

between a fuse and a breaker : 

j 

BRKBEK - calculates coordinat.^on indices ; 

between two breakers 

RELFUS - calculates coordination indices ? 

ji 

between a relay and a fuse | 

RELBRK - calculates coordination indices j 

between a relay and a breaker i 

The above subroutines calculate the following 

i 

probabilistic coordination indices Al, A2 and RSC, which ; 

are defined in the subsequent section. They use the set 
of subroutines (a) to characterize the protective elements 

5 

and the Short Circuit Program to determine available I 

fault currents* 

(c) PROTECTION SUBROUTINES ^ 

Two subroutines have been developed to check 

protection: j 

CHEGKH - checks if the system is 

protected against overheating 

- checks if the system, is pro- ' 

tected against motor starting currents. 


CHECKS 


These require data on heating characteristics of network 
elements and motor starting current vs time requirements. 

The above set of subroutines uses the set of subroutines (a), 

(d) PLOTTING SUBROUTINE 

This division includes a plotting program consisting 
of two parts: 

a. the routine PLOT which creates the characteristics 
to be plotted in digital form 

b. the part which utilizes the graphics system 
software to create the graph (Datagraphics 4460 ). 

A modification of the above routine under the name subroutine 

PLOT is very useful in connection wiuh the optimizing program. 

Note: This set of subroutines utilizes the device 

characteristic subroutines and the short 
circuit program. 

(e) COORDINATICM OPTIMIZATION SUBROUTINES 

.These subroutines provide a computerized optimization 
procedure. The designer has to determine all the candidate 
settings which can be assumed by the system. The program 
calculates an overall coordination index for every combination 
of settings which yields the higher overall coordination 
index. The format of the required data is explained under 
the program FBBBOG, 

The program FBBBOG refers to a system which includes a 
fuse and three breakers. This is the case of a typical 
substation at the Cape Kennedy Power System. 


The program can be used as the basis for other programs 
which treat other cases, for example a substation implemented 
by a primary fuse and two breakers on the secondary. Actually 
this modification is offered as an option of the program FBBBOC. 
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2. DEVEiiOPMENT OF THE GOOEDimTlON PROCEDTJEE 

The following pages present the approach that has been taken to the 
coordination problem. This approach provides two main features: 

(1) An algorithm which determines the optimum coordinated settings. 

The algorithm is employed in a computer program which determines 
the optimum coordinated settings 

(2) A plotting routine which will graph on a common set of axes 
the characteristic curves of the protective devices employed 
by the system. 

The optimizing algorithm applies to a radial distribution system 
such as the one at KSC. Looped systems are not considered. 

2.1 REPRESENTATION OF DEVICE TIME VS. CURRENT 
CHARACTERISTICS 

The first problem encountered in system protection coordination 
is how to reproduce the operating characteristics of the protective 
devices. Basically an approximation of the function 

t = £(I) 

is required. 

To achieve this, it has proven practical to perform a coordiviate 
transformation to a logarithmic scale as follows: 

t' = Cj^ log t + 

X = log I + 
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where and constants serving normalization in such a way that 

0 ^ 1 
0 ^ X s; 1 

In this case, the function t == f(I) to be approximated becomes 

t ^ - g (x) , (0 :S X s: 1) , (0 :£ t ^ s; 1) 

The scheme used to approximate the function 

t* - g(x), (0 s£ X £ 1), (0 s: t' 5 1) 

is a piece-wise linear approximation of the curve; that Is, the curve 
t’ = g(x) is approximated with pieces of straight line. This scheme is 
very •'-'mple and gives good accuracy. The accuracy depends on the number 
of line segments chosen; increasing as this number increases. For a 
typical breaker trip device curve, thirteen points give accuracy greater 
than 3%. 

The following figure illustrates the method used. It should be 
clarified that in practice the points are chosen in a very dense mode 
where the curve is highly curved and in a dilute mode where the curve is 
close to a straight line. Some specifications define maximum clearing 
time curve and minimum clearing time curve. In tliis case, each curve is 
treated separately. 

This approximation scheme provides the basis for the developed sub- 
routines which are given in Appendix 1, Each subroutine can provide the 
maximum clearing time and the minimum tripping time of the corresponding 
protective device if the characteristics (settings or size) of the device 
and the current are determined. Subroutines are in.iluded for almost all 
devices employed in the KSC network. 



The above subroutines are used in the plotting program to 
graphically represent the characteristics of the devices which operate 
in a system under study. The characteristics are referred to a common > 
current axis and then plotted. Details are explained in the plot routine. 

next probabilistic concepts of coordination which will he proven 
to be the analyf;ic tool for an optimizing procedure are introduced. 

■ 2.2 OPERA.TIONAL CHAEACTERISTICS OF DEVICES 
In connection with the time-current characteristics of the devices, 
some operational factors are of importance in attempting to provide 
computer assistance to coordination. These pertain mainly to tolerance 
values on the specified operational curves. The specified time to trip 
for a certain level current is not a fixed value and varies due to a 
variety of reasons. This suggests that the time to trip of a device 
under sustained and fixed current level is distributed with some pro- 
bability around the manufacturer's specified value. There are several 
conti^^utions to this distribution. The following items Include the 
most important effects: 

1. Overtravel (applies to induction type relays). 

2, The burden of the protective device circuit has an error of 
its own which introduces uncertainty about the value of the 
operating quantity (current or voltage). 

3, The dead time (the time to clear the arc in the breaker) is 
not constant. Its duration depends on the condition of the 
breaker, the relative time to the cycle in which the arc 
starts, and breaker interrupting rating. 

4. The fault current waveforms are usually distorted and therefore, 
the conditions of operation differ from those used for plotting 
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the time current curvB. 

5. Ambient temperature has a great influence especially on thermal 
trip unit devices. In this case, the thermal history of the 
system before the fault occurrence is significant. 

6. Mechanical tolerances introduce an error especially on devices 
with movable parts and on filament type fuses, 

7. Operator Introduced error (important for devices with continuous 
settings). 

In addition, the computer approximate representation scheme 
introduces an error. 

Most of these sources of error may be asstraied to have Gaussian 
distribution. Assuming that each distribution has a standard deviation 
then the total error will be. Gaussian distributed with a standard 
deviation: 



with summation over all applicable i. 

The assumption of the Gaussian distribution is made witn reference 
to the t* variable (t* “ t + C* ) which is in agreement with the 

mathematical requirement for the Gauss distribution to have a domain 
_« <t' <«, Note that the model is physically consistent since the time 
axis o<t<« transformed to the t' axis becomes -03<t'<«. 

The values of which apply for a given device are not all 
precisely known. Furthermore, the computational details connected with 
the Gaussian distribution are involved. Therefore it is practical to 
assume a less computationally involved distribution. As such the uniform 
probability distribution has been accepted. This distribution may be 
readily defined since the maximum and minimum tripping times are generally 
known. Using these the error distribution is determined as shown below. 
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2.3 MEAS’jRES OF COORDINATION 

In order to utilize the computer to assess coordination between 
devices, it is necessary to define a quantitative measure of coordination 
to be used as a basis for comparison. These measures should provide in 
some sense an indication of the difference between interrupting times at 
a given current level of the devices being compared. As discussed pre- 
viously, interruption time is regarded as a random number distributed in 
a uniform fashion about a nominal value. It therefore makes sense to 
develop a statistical measure to define coordination. In the following, 
the concept of coordination reliability is introduced. This is then 
extended to develop one additional index of coordination. 

A, Reliability of Coordination 

Consider a I'adial circuit with two protective elements M and DB. 
Between them is an arbitrary network section which may or may not include 
transfomers. In general, if I is the current flowing through device 
DB, then 1^^ = current flovzing in the device DA, Define 

r as the transformation ratio. 
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Suppose a fault current I through this radial is sustained until 

one of the devices trips. To assess coordination between DA. and DB, it 

is necessary to know a priori which element interrupts first. 

Assume that device specifications indicate that under sustained 

current I the interruption times for DA and DB are ^*Amin 

device DA. and t’ » t’ . for device DB, According to the previous 
Bmax Bnn.n 

assumption of uniform distribution, the interruption times of the devices 
will then be distributed as : 


probability] 
density 
function 




^ A I? _i. 


"Bmin AnjfLn Bmax 

i > 

I ' 

V' 


' 

Device DB will interrupt in the interval (x^, + dx^) with 

probability 


^Amax 


dP„ ^ 


B t' - t' . 

• Bmax Bmin 


dx. 


B * 


Device DB will act first if device DA interrupts in any time x 


where x>x . The probability of this event is; 
B 




t' - t' * t’ “ t' 

Bmax Bmin Amax Amin 


if x„ > t* 


t* - t’ 

Bmax Bmin 


^ Amin 

dKs if 




(t*. " X ) 

Amax B 


B. Expected Reliability of Coordination 


Another important measure of system coordination relates to 
coordination on the average over the range o£ possible fault current 
levels . 

The range of possible fault currents of a protective device is 
defined by the c?nrrents in the interval from the pickup value up to the 
short cicruit capacity at the point of application of the protective 
device . 

Now consider the following situation. Along the radial there are 
two adjacent protective devices DA and DB, and between them there is an 
arbitrary network section. 



At the position of the fault, the rated current is and the short 


DB 


circuit capacity is 


.DB 


The device DB is set to trip if the current is greater than 

DA 

Similarly, for device DA we have 

It is obvious that the common range of possible fault current for 
the two devices is 

■^Pickup Short 

For every current we can calculate, the reliability of coordina 

DB 

tion between the devices DA and DB, R 
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Therefore, the probability that device DB will interrupt first is: 

tV 


Bmax 


dp. 


R - r 

^ Bmin 

The probability number R is called the reliability of coordination 
between devices M and DB for the fault current I. 

This probability is a function of the current, and the specific 
settings of the devices. 


R = R (I, S^, S^, r) 

where S. CS_) represents the settings of the device DA. (DB), and r is the 

A. B ' 

transfoimation ratio. 

It can be shown that:. 

R (I, S^, r) = 1 - R (I, S^, S^, r). 

It should be noted that this figure gives a measure of reliability 
at only a single level of fault current. One particular value of fault 
current for which the reliability index is significant is at the short 
circuit capacity of the associated line. This is due to the fact that 
for an underground cable system such as the system under consideration, 
most faults are expected to draw fault current near the short circuit 
capacity. The reliability index for this point is: 


This number then indicates coordination between elements DA and DB 
at one particular current value and is of great importance. 
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Then, on the average we will have 
E[R (Ids, Sb> r)] 

where the expectation extends over the conraon range of fault currents. 

This measure of coordination is designated as A2. If the distribu- 
tion of fault currents in the range of lj^g» f known, then A2 

may be computed explicitly as 


all I. 


DB 


The probability density distribution f(I^g) is unfortunately not 
readily determined. 

It: has been determined, however, that in situations where the 
coordination is reasonably good, the measure A2 is not very sensitive 
to the specific distribution. It is therefore reasonable to make the 
simplifying assumption that fault currents in the range C^j^ckup*^Short^ 
occur equally likely. That is 


L 




jDB __ ^DB 
Short Pickup * 


.TB 


Pickup ^ ^ 


DB 

DB " "Short 


otherv/ise 




Then 


A2 = 


.1 


DB ^DB 

Short Pickup 


.DB 

short 

^DB 

Pickup 


^ ^^DB’ ^A’ 



h’ 
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Interpret: at Ion of A2 as a distance measure 
The previoxisly defined coordination measure A2 can also be 
interpreted as a measure of the relative position of the two cutrves. 




t' = t' 

B B 




For example, if A2 is very close to unity we can conclude that curve 

above the cuve t'^ = 

Guided by this interpretation a more rigorous measure of coordina- 
tion may be developed. Suppose the minimum clearing time curve of device 
DA and the maximum clearing time curve of device DB are 


raxn mun 


and t'^ = respectively. 


max max 

Another measure of coordination is then defined by the number A1 


computed as 


A 1 = E[RC(f^ - f 3 (I33)), S , S3, r)] 

mxn max 

where the expectation is taken in the same sense as for A2, 

Obviously, if number A1 has a value close to one, the minimum 
clearing time curve of device DA is well above the maximum clearing time 
curve of device DB and therefore that coordination exists between the 
two devices. 

Thus two different coordination measures are defined. The number 
A2 represents a physical quantity: It gives the probability of coordinated 

operation of the two devices over the entire operating current range. 

Its value provides information about possible overlaps of the associated 
time delay bands. 

The number A1 describes the relative position between two curves. 

For example, if A1 = 1.0, it is concluded that the two curves never 


intersect and the distance between them is everywhere greater than a 
prespecified value. Therefore, A1 is a stricter measure of coordination 
and can be considered as a safety factor. 

2.4 COMPUTERIZED GALGUIATI0I5 OF THE 
COORDINATION INDICES 

The subroutines RELFUS, RELBEK, FUSBRK and BEKBEK given in 
Appendix 2 calculate the coordination indices A2, A1 and R for the 

SC 

pairs of devices Relay-fuse, Relay-breaker, fuse-breaker and breaker- 
breaker respectively. 

Each subroutine originally defines the set of all possible fault 
currents as it is described in the particular subroutine. Then it 
divides this set into NN intervals. It calculates the reliability of 
coordination in every part. Finally the coordination index is the sum 
of the calculated reliabilities weighted with the length of the corres- 
ponding interval. 

The described procedure is applied for the indices A1 and A2. 

Last the reliability of coordination at the short circuit capacity RSC, 
is calculated. 

2.4,1 Stunmary 

In the preceeding text three coordination indices have been defined, 
A2, Al, and RSC. Tbeir physical meaning is as follows: 

A2: It assumes values between zero and one and represents an 

expectation of good coordination between two devices over a 
predefined range of fault currents. If its value is 1.0 it is 
to be concluded that the coordination is very good over the 
mentioned range of fault currents. 
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It is obvious that in a system where any kind of fault 
can occur with any value of fault impedance, the index A2 
is very important. This is so, because under these conditions 
the level of the fault current can be anywhere above the load 
current of the system and below the bolted three phase or line 
to ground fault whichever is greater- The above statement should 
be modified accordingly if considerable asyinmetry of fault 
currents is present. In other words, values of A2 near unity 
guarantee selectivity over the range of classified fault 
currents . 

Al : It assumes values between 0 and 1.0, and represents a safety 
factor of good coordination between two devices over a pre- 
defined range of fault currents. 

Al expresses an average of safety factor over the mentioned 
range of fault currents. Therefore, values of Al close to 
unity guarantee selectivity with a safety factor not only at 
- one level of fault current but over *j. classified fault currents. 

A comparison between A2 and Al reveals that Al is a stricter 
coordination index in the sense that if Al = 1.0 then iEunediately 
A2 * 1.0, But if A2 = 1.0 then is not sure that Al=1.0. 

Al is a very useful index in characterizing the reli- 
ability of a system’s selectivity. Or when xincertainty for 
the trip devices has been acctimulated by the pass of time due 
to the lack of field tests. 

ESC; It assumes values between 0.0 and 1.0 and expresses the reliability 
of coordination between two devices at the short circuit capacity 


of the system. 
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Unlike A1 and A2 it provides information at one and only 
one level of fault currents. The usefulness of this index is 
due to the fact that for some systems (e.g. underground cable 
or aerial cable) the levels of fault currents at a point of the 
system are confined with high probabi lity in a rather narrow 
area just below the short circuit capacity of the system at 
that point. Given that the time vs current characteristics 
of the protective devices are smooth over wide intervals it can 
be concluded the RSC represents the coordination status of such 
systems with high probability. 

Fom the preceeding discussion it is obvious that each eoct dination 
index can be used for appropriately selected purpose of coordination, for 
example RSC can be used to assess coordination at the short circuit 
capacity of the system. 

In general, for comparison purposes, it is expedient to express 
coordination with one number which will reflect all three defined co- 
ordination indices. The linear combination of Al, A2 and RSC has been 
adopted as follows: 

RELOC = xA * A2 + xB * Al + xG * RSC 
where xA + xB + xG = 1.0 
and therefore 0.0 s: RELOC ^1.0 . 

The coefficients xA, xB, and xC provide the flexibility to the 
designer to accomplish his own tendencies to the solution of the coordina- 
tion problem by appropriately selecting the values of xA, xE, and xC. 

For example, suppose that designer A desires to have the protective 
devices coordinated at the short circuit capacity of the system. Then it 
is sufficient to select xA = 0, xB = 0, and xC = 1.0, 


Now the number RELOC can be used as the comparison criterion between 


two different solutions to any coordination problem. 

The above idea leads to the optimization algorithm. This algorithm 
will be presented for two devices: 

Step 1: All the candidates settings of the devices/ or 

devices are defined by the designer. The device which is supposed 
to trip first is specified. The short circuit capacity of the 
system at the point of application of the other device is 
calculated. Then reasonable range of possible fault currents 
is defined. 

Step 2 


Step 3 


Step 4 


Step 5 


Step 6 


Step 7 


Step 8 


Assume a solution i.e. devices and/or settings 
Compute Al, A2, and RSC for tha above solution 
Compute EELOC = xA * A2 + xB * Al + xC * RSC 
Register above solution as "optimum" 

Assume a different solution 
Compute Alj A2, and RSC for this solution 
Compute EELOC” xA * A2 + xB * Al + xC * RSC 
Step 9: Compare the two EELOC values. If second value is 

greater, register the current solution as "optimum." If 
first value is greater do not change the "optimum" solution. 
Finally keep the higher value of EELOC as comparison 
reference for the next comparison in this step. 

Step 10: Asstnne another solution, different than the previous. 
Then go to Step , Incase that all the candidate solution 
have been used, terminate the loop and go to the next step. 

Step 11: Print out the currently registered "optimum" solution 

and relevant information. 
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The optimization program expands the above algorithm to more 
than two devices. The selected values of the coefficients xA., xB, and 
xC are: xA. = .35, xB = .15, and xC == .50, 


2,3 COMPUTERIZED PLOTTING 

The development of the device characteristic subroutines makes 
easy the plotting of different characteristics of protective devices by 
the computer. The procedure is very simple. Suppose there are two 
devices A and B on a radial. Also, suppose a fault current flows 
through device B. The time to operate t„ can be calculated by calling 
the appropriate subroutine of this device. Now when fault current 
flows through device B, fault current flows through device A. can 
be computed with the current transformation ratio between these two 
devices. The time to operate t, , of device A can be determined by calling 
the appropriate subroutine for this device. Now on the predefined common 
coordinate system, say t versus the computer puts the points Ig) 

and (t^, Ig). The above procedure is repeated for many values of 
over the range of possible fault currents Tg and finally a line is drawn 
through all the points (t^, 1^) and another through all the points (t^, 1^), 
Note that t. (or t_) can represent minimum operating time or maximum or 
their arithmetic mean, 

A program has been developed which expands the described procedure 
tc more than two devices. 


Ov'erlayed characteristics on a common coordinate system of the 
devices in a system is the best information about the coordination status 
of the devices. Therefore the plotting program can be used as a check 
to the output of the optimizing program (in fact the optimizing program 

. y . . ' 
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plots the optimma solution) or to evaluate by inspection any manually 
given solution to the coordination problem. 

3, APELICA.TIONS 

In this section, application of the coordination and plotting 
programs are covered in detail, A. specific example of application to a 
typical KSC substation configuration is examined initially. This example 
will illustrate the procedure for most existing KSC installations and a 
particular subroutine has been developed to treat this configuration 
(double fed substation with primary fuse, main secondary breaker, tie 
line breaker and feeder breakers) . Following the example procedures 
necessary to accomodate alternative substation configurations are 
explained. 

t 

As an example, the coordination of the substation 804 located in 
the VAB is studied. This substation is a typical double fed substation 
in the system. The next figure shows the substation and its interconnec- 
tion with the rest of the system. 

The study is worked out in two steps! 

a. Preparation of the data cards for the plotting program is 
carried out for the devices and settings as they exist. This 
results in a computer plot of the overlayed protective device 
characteristics. From the output the user can make an evalua- 
tion of the coordination status by inspection. 

b. Preparation of the data cards for the optimizing program is 
carried out with comments which will be helpful for other cases. 

The optimizing program also utilizes the plotting to provide 
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a plot of the overlayed time vs current curves of the devices 

for the optimal settings. 

3.1 USING THE PLOTTING ROUTIKS 

Use of the plotting routine is exemplified using substation 804 in 
the VAB. Data card preparation and output format are discussed. 

The program has been designed to accomodate at most seven devices, 

2 relays, 2 fuses and 3 breakers. The following identification names 
have been assigned to them: relay A, relay B, fuse A, fuse B, breaker A, 

breaker B, and breaker C, 

The program requires eight data cards. The first one should con- 
tain general information about the system. The remaining seven correspond 
one to one to the seven devices. A missing device is identified by a 
blank card. 

Preparation of data has been made under the following assumptions; 

(1) Only fault currents are considered. The load current of 
interconnected lines and apparatus is ignored. Current trans- 
formation ratios are computed under this assumption 

(2) The interrupting capacity of the breakers is assumed to be 
20,000 Amperes, This is due to lack of information and 
should be modified accordingly when this information is 
available, 

(3) Current transformation ratios that cannot defined (for 
example between a non-existing device and another device) 
should be given the value 1.0, 

It is desired to plot the characteristic curves of the protective 
devices in substation 804 located in the VAB. 
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fuses; type EJ 

2 ® ® 


750 kVA 

A/y^ 

"B” 


© © P 


A/W 

750 kVA 

t»A» 



VAB 
SS 804 






breakers ; type AK by GE 


This requires that the program plot the characteristic curves of 
the fuse (on the primary), of the main breaker, the tie line breaker and 
any one of the feeder breakers. 

Therefore, data for one fuse and three breakers is required. The 
data cards for the relays A & B and fuse A should be blank. 

Preceding all data cards the user should provide one card on which 

appears the title to appear on the plotted output.' 

1st data card 

1-10 Current transformation ratio between 

relays A and B. (Relays A and B do not 
exist, therefore insert 1.0) 

11-20 Current transformation ratio between 

relay B and fuse A. (Relay B and fuse A 
do not exist, therefore insert 1.0) 

21-30 Current transformation ratio between 

fuse A and fuse B. (Fuse A does not exist, 
therefore insert 1.0) 

31-40 Current transformation ratio between fuse 

and main breaker is 13,800/480 = 28.65 

41-50 Current transformation ratio betvjeen main 

breaker; insert 1.0 and tie line (worst Case) 

51-60 Again insert (1.0) (worst case) 



26 


% 


O' 


61 Insert 1 since substation 804 belongs to 

the Launch Complex 39 (Launch Complex 39 
is record 1 on the network data tape 
used by the KSC short circuit program. ) 

2nd data card 

blank (relay A does not exist) 

3rd data card 

blank (relay B does not exist) 

4th data card 

blank (fuse A does not exist) 

5th data card 

1-5 insert the symbolic name of the bus where 

the fuse is located, i.e. A24 

6-10 insert in column 10 the integer 1 since 

the fuse is of type EJO-1 

11-15 the code for the size of the fuse. Size 

of the fuse 80E. From subroutine FUSEEJ 
the code for this size is 15. Therefore 
insert in columns 14 and 15 the number 15 

16-25 the current rating of the fuse, i.e. 160,0 

6th data card 

1-5 insert the symbolic name of the bus where 

the main breaker is located; that is C41 

6-15 the rated pickup current of the breaker, 

that is 1200,0 

16-25 the interrupting capacity of the breaker, 

(Assumed to be 20000.00 anips.) 


27 


26-30 the code for the type of the breaker, 

that is 1 (breaker type AK) 

To fill in the rest of this data card, reference must be 
made to subroutine BiEiKRAK 

31-35 the code for the series trip device. 

Insert 1 (EC-1 series trip device) 

36-40 code for ube long time delay band. 

Max — )■ 3 

41-45 code for the short time delay band. 

Intermediate — j- 2 

46-55 long time coil rating in per unit. In 

our case equals 1.20 

56-65 short time pickup setting as a multiple 

of the rated pickup value. In our 
case 5.0 

7th data card 

The format is imilar to the 6th card. Insert the numbers 
corresponding to the tie line breaker as follows: 


1-5 

C4l 

6-15 

1200.0 

16-25 

20000.0 


26-30 


1 



46-55 


1 

3 

2 

0.8 


56-65 


2.5 
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8th data card 

Select any one of the feeder breakers since each one is 


set in an 

identical way. 

1-5 

G41 

6-15 

300-0 

16-25 

20000.0 

26-30 

1 

31-35 

1 

36-40 

2 

41-45 

1 

46-55 

1.0 

56-65 

6.0 


Output 


i'he curves for the seven devices are 

plotted 

on a comnion scale. The 

lines are distinguished by the length and width of 

the dashes and by the 

nuiaeric character 

printed on 

each point as 

follows: 


Device 

Uumber 

Dash Length 

(rasters) 

Line Width (rasters) 

Relay A 

1 

solid 


4 

Relay B 

2 

solid 


8 ■ 

Fuse A (Upper) 

3 

32 


2 

(Lower) 

3 

32 


4 

Fuse B (Upper) 

4 

32 


8 

(Lower) 

4 

32 


16 

Breaker A (Upper) 

5 

64 


2 

(Lower) 

5 

64 


4 

Breaker B (Upper) 

6 

64 


8 

(Lower) 

6 

64 


16 

Breaker C (Upper) 

7 

256 


2 

(Lower) 

7 

256 


4 
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3.2 USING THE OPTIMIZING ROUTINE 
In this section the optimizing procedure is employed on Substation 
804 in the VAB. The example outlines data preparation and procedures for 
applying the routine. The alternative substation configurations are 
given. 

General Observation 

The closest protective device to the substation is the breaker at 
the SWGR 700 (or SWGR 702 at the other side). Consideration is given to 
one side only (SWGR 700) since the other srde is similar. The circuit 2 
of the SWGR 700 which feeds transformer A in substation 804 is protected 
with an oil circuit breaker governed by an lAC relay. This circuit 
feeds a number of substations and therefore the relay's pickup current 
is well above the rated current of the transformer A in substation 804. 
The pickup current of the relay (circuit 2 at SWGR 700) is 480 amperes 
and the rated current of transformer A in substation 804 is 31.4 amperes. 
Under these conditions it is certain the coordination exists between 
the relay and any protective device in substation 804 (under the logical 

assumption that only one fault occurs in the system at a time). 

and the rated current of transformer A in substation 804 is 31.4 amperes. 

Under these conditions it is certain that coordination exj sts between 
the relay and any protective device in substation 804. 

Therefore it is necessary only to study coordination of the pro- 
tective devices which are located in. the substation and after it. In the 
case of the substation 804 it is required to determine coordination of 
the primary fuse, the main secondary breaker, the tie line breaker and 


the feeder breaker. 
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The Mathematical Model 

The figure shows the substation and the symbolic names of the 
buses according to the model (short circuit study of the system) 



is located. The other side is similar. 

1. In the analysis presented here, the worst case of coordination is 
considered: 

Suppose transformer B is out. Then, in order not to interrupt 
service at the load of the transformer B the tie line breaker has to close 
and the main secondary breaker of transformer B to open. In this case 
the following radial system exists: primaiy fuse - transformer '*A” - 

main secondary breaker - tie line breaker - anyone of the feeder breakers 
on side *'B". 

Preparation of the data has been made under another assumption, 
i.e. load currents have been neglected. Therefore the current transforma- 
tion ratios represent fault current transformation ratios. 
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Data Preparation 

Step 1; Preparation of the data cards 

Preceding all data cards the user should provide three cards in 
which he can punch anything he wishes to be the title of the program's 
output . 


In this case a good title can be; 

First card: SUBSTATION 804 LOCATED IN V.A.B 

Second card: TEAKSFOEMER A SIDE 
Third card : Blank 


1. first card 


Position 


1-6 

7-12 

13-18 

19-24 

25 


Insert 

A24 

G41 

C41 

C41 


Explanation 

bus name at the location of the fuse 

bus name at the location of the main 
secondary breaker 

bus name at the location of the tie 
line breaker 

bus name at the location of the 
feeder breaker 

since we do not have data on the over- 
heating of the transformer (at the 
present time) 

since the substation does not feed 
any motors. 
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Position Insert 

27 1 


28 


0 


Explanation 

The «ubstion belongs to the Launch 
Complex 39. (Launch Complex 39 is 
record 1 on the network data tape 
used by the KSC short circuit program) 

Double fed substation with primary 
fuse, main secondary breaker, tie 
line breaker and feeder breakers. 


2. Second card 

The interrupting capacities of the breakers are required. 
At the present time this information is not available. 
Therefore the value 20,000.00 in amperes is assumed: 


Position 

Insert 

1-10 

20000.0 

11-20 

20000.0 

21-30 

20000.0 


3 . Third card 
Position 
1-10 

11-20 

21-30 

31-40 


Insert 


28.65 


1.0 


1.0 


1200.0 

1200.0 


Explanation 


13800 

480 


= current transformation 


between fuse and main breaker 


current transformation ratio between 
main breaker and tie line breaker. 

current transformation ratio between 
tie line breaker and feeder breaker. 

rated pickup current in amperes of 
main secondary breaker 

rated pickup current in amperes of 
tie line breaker. 


41-50 


Position 


Ins ert 


Explanation 

51-60 300.0 rated pickup current in amperes of 

feeder breaker. (Rote: We usually 

choose the feeder with the higher 
normal load. In our case every 
feeder is designed for the same 
normal operation load. ) 

Fourth card 

This card refers to the primary fuse which is of type EJO-1. 

Therefore subroutine FTJSEEJ should be referenced 


Position 

Insert 

Explanation 

1-5 

14 

fuse type EJO-1, 65E, size DD 

6-10 

16 

fuse type EJO-1, lOOE, size DD 

11-15 

1 

step of the do loop. 


At the present time the primary fuse at the substation 804 
is of type EJO-1, 80E, size DD. Looking at the corresponding 
subroutine (FUSEEJ) it is seen that three fuses of size DD are 
utilized, namely 65E, 80E, lOOE with the codes 14, 15, 16 
correspondingly. These three fuses can be interchanged without 
any change in the system. In order to nominate all three fuses 
as candidates insert the number 14, 16 and 1 as has been described. 

Suppose it is desired to retain the fuse 80E in the system. 
Then insert 15 in the two positions; 1-5 and 6-10 and leave 
everything else as it is. 

Position Insert Explanation 


16-20 


1 


the code for the subroutine FTJSEEJ 
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Fifth card 

This card refers to the main secondary breaker which is of 
type AK. Therefore subroutine BjRKEAK should be referenced. 


Position 


11-15 


16-20 


21-25 


26-30 


31-35 


36-40 


41-45 


Insert Explanation 

1 the main secondary breaker is of type 

^ AK with EC-1 series trip device. 

Looking at the subroutine BKKRAK it 

1 

is seen that the code for the EC-1 
series trip device is 1, Furthermore 
it is not desired to change the 
series trip device. So the program 
is instructed to consider the series 
trip device with code 1 that is EC-1. 

1 code for the long time delay band. 

2 The program will consider the minimum 
band (code 1), the intermediate band 

^ (code 2) and the maximum band (code 3), 

For the interpretation of the code 
see subroutine BRKSAK 

1 code for the short time delay band. 

2 Same comments as for the long time 
delay band can apply. 


46-50 


51-55 


the long time coil rating expressed 
in percent. We want the computer to 
consider the settings between 80X 
and 120% with step of 20% that is 
the settings 80%, 100% and 1207o. 


56-60 


20 


35 


Position 

61-65 

66-70 

71-75 


76-80 


Insert 

3 

5 

1 


1 


Explanation 

short time pickup current. From 
the manufacturer data it can be seen 
that the short time pickup current 
can take any value between 2 and 5 
times the rated pickup current of 
the unit. Here the computer is to 
consider all the values between 3 
and 5 times the rated pickup current 
(i.e. 1200.0 amperes) proceeding 
with step 1. In other words the 
following values 3, 4 and 5 times 
the rated pickup current - 

code for the breaker type AK. The 
code refers to the subroutine BRKRAK. 


6. Sixth card 
Position 
1-5 
6-10 
11-15 


Insert 

1 

1 

1 


Explanation 


EG-1 series trip device 


16-20 1 

21-25 3 

26-30 1 

31-35 1 

36-40 3 


minimum-intermediate and maximum 
long time delay band. 


41-45 


1 


minimum- intermediate and maximum 
short time delay band. 
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Position 

Insert 

Explanation 

46-50 

80 


51-55 

120 

80%, 100% and 120% long time coil 
rating 

56-60 

20 


61-65 

2 

2, 3 and 4 (times the rated pickup 

66-70 

4 

current) short time pickup settings 

71-75 

1 


76-80 

1 

breaker type AK and therefore 
subroutine BRKRAK. 

Seventh card 

This card refers to the 

feeder breaker No. ICB4 - lOCLA which 

is an AK-2A-25 

breaker. 

Therefore the subroutine BEKRAK should 

be referenced. 

Preparation of data and explanation is similar 

as per the fifth data card. 

Position 

Insert 

Explanation 

1-5 

1 


6-10 

1 

EG-1 series trip device 

11-15 

16-20 

1 

rainimmn, intermediate and maximum 

21-25 

3 

long time delay band 

26-30 

1 


31-35 

1 

minimum, intermediate and maximum 

36-40 

3 

short time delay band. 

41-45 

1 
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Position Insert 


Explanation 


46-50 

51-55 


56-60 


100 

We desire to keep the rating of the 
long time coil fixed at 100 percent, 

20 


61-65 

66-70 

71-75 


76-80 


4 from the manufacturers data it is 

seen that the short time pickup 

8 

current can be any number between 4 

-r 

and 10 times the rated pickup current. 
Therefore the computer is to consider 
only the settings 4, 5, 6, 7 and 8. 

1 type AK breaker. 


8 . Eighth card 

Position 

1-4 


Insert 

PLOT 


Explanation 

The word ”PL0$" here causes the 
optimal coordination curves to be 
plotted by the computer. This 
card is omitted if no plot is 
desired. 


MOTES 

1. When an integer number is to be puiiched in a certain 
position bear in mind that the computer will interpret a blank 
space as zero. Therefore when it is required to punch the 
integer 14 in the columns 1 through 5, al should be punched 

in column 4 and 4 in column 5, 

2, The program can handle as well single fed substations which 
are implemented with primary fuse, main secondary breaker and 
feeder breakers. This is offered as an option. 
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ImplemeTitatiorL of the above modification (option) is achieved 

by punching a one in column 28 of the first data card. In this 

case the computer skips the part which refers to the second 

breaker and the optimiaing algorithm applies to the fuse and two 

breakers only. For practical purposes (to avoid warnings) the 

data card which refers to the second breaker should be identical 
to the data card which refers to the first breaker. Finally 

any output referring to the second breaker (breaker B) should 
be ignored. 


PART II. Appendices. 
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INDEX OF SUBROUTINES 

Appendix I. 

Sygibollc name 
RELIAC 

RELACO 

BRKRAK 

BRKELA 

BRKROD 

BRKRDB 

MCCBGE 

FUSEEJ 

FUSCLE 

FUSECL 

FUSESM 


Description Page 

simulates the characteristics of 
relays type lAG by G.E. 

simulates the characteristics of 
relays type CO by Westinghouse 

simulates the characteristics of 
breakers type AK by G.E, 

simulates the characteristics of 
breakers type LA by Allis- Chalmers 

simulates the characteristics of 
breakers type OD by I.T.E. 

simulates the characteristics of 
breakers type DB by Westinghouse 

simulates the characteristics of 
molded case circuit breakers 
by G.E. 

simulates the curves of fuses type 
EJ-1, EJO-1 and EJ-2 by G.E. 

simulates the curves for fuses type 
CLE-1 and CLE-2 by Westinghouse 

simulates the curves of fuses type 
CL-13 by I.T.E, 

simulates the curves of fuses type 
SM by S 6: C Electric Company 

simulates a piecewise linear func- 
tion defined by K points C2dN^0), 
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Appendix 2 


Symbolic name 
FUSBRK 


BRKBKR 


RELBRK 


RELFUS 


Description Paae 

computes the coordination indices 
Alj A2 and RSC between any fuse 
and any breaker existing in 
Appendix 1. 

computes the coorindation indices 
Al, A2 and RSC between any two 
breakers existing in Appendix 1. 

computes the coordination indices 
Al, A.2 and RSC between any relay 
and any breaker existing in 
Appendix 1. 

computes the coordination indices 
Al, A2 and RSC between any relay 
and any fuse existing in 
Appendix 1. 


Appendix 3 

Symbolic name 
CHECKH 


CHECHS 


Description Page 

it checks if a prespecified fuse 
(from the existing in Appendix 1) 
clears any fault before any 
damage can happen to the system 
due to heat dissipation 

it checks if a prespecified 
breaker (from the existing in 
Appendix 1) does not operate 
when an inrush motor starting 
current occurs in the system. 
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Appendix 4 

Main program 
PLOT' 


Subroutine 

PLOT 


Appendix 5 

main FBBBOC 


Page 

it overlays on a common coordinate 
axis the time vs current charac- 
teristics of at most seven devices 
(2 relays, 2 fuses and 3 breakers) 
from the inventory in Appendix 1, 

it overlays on a common coordinate 
axis the time vs current charac- 
teristics of a fuse and three 
breakers frcsn the inventory in 
Appendix 1. 


it provides the optimum coordinated 
settings of the protective devices 
in a substation implemented by a 
fuse, a main secondary breaker, 
a tie line breaker and feeder 
breakers . 

An option of this program 
handles a substation implemented 
by a fuse, a main secondary 
breaker and feeder breakers. 


APPENDIX X 


DEVICE CHAEACTERISTICS SUBROUTINES 


These subroutines simulate the time vs current characteristics of the 
protective devices employed in the Kennedy Space Center system. 

The piecewise linear approximation of a curve is generally used. Each 
curve is defined with a number of points which differs from case to case. 
The appropriate characteristic is determined with minimum device 
paramet^its. 

The subroutines can be used in two ways; 

A. To return the operating time of the device at a given current, and 

B. To retrieve the time vs current characteristic of the device 
for predetermined settings. 


SUBROUTINE RELIAC 
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SUBROUTINE HELIA.C 
General Descri-ptxon 

RELIAC is a FORTRAN subroutine which returns the time required by 
a relay type lAC of General Electric to close a tripping circuit when 
the settings of the relay and the current are specified. 


Calling Sequence 

CALL RELIAC (XR,NTCHR,TDS , YR) 

where 


XR is the value of current expressed in multiples of the pickup 

value. 

NTCHR is a coda for the type of the time characteristic curve and 
assumes integer values as follows: 

NTCHR =1 = inverse time characteristic 

2 = very inverse time characteristic 

3 = extremely inverse time characteristic 

4 = inverse, short-time characteristic 

5 = inverse, long- time characteristic 


TPS is the time dial settings and assumes real values in the 
field (.5 - 10.) 

YR is the variable which will assume the value of the time to 
close the tripping circuit in the transformed system of coordinates. 
Therefore, the same time in seconds will be; 

time to close the tripping circuit = .01*(10.0**(5.0*YR)) 


t 
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Coimnents 

For each family of curves (for example very Inverse time character- 
istic), the curves corresponding to the time dial settings .5, 4.0, 
and 10.0 ■^lave been stored (with 13 points each) under the symbolic names 
and (XC,YC) respectively. 

Analytically; 

CXC(K),YA(I)) , IDS - .5*^1 

(XC(K),YB(I)) , TDS = 4. i K = 1,13 

(XC(K),YC(I)), TDS = 10. ] 

I = 1,13 = inverse time characteristic 

I = 14,26 = very inverse time characteristic 
I = 27,39 = extremely inverse time characteristic 
I 40,52 = inverse, short-time characteristic 
I = 53,65 = inverse, long-time characteristic 

Curves for other time dial setting are approximated by interpolation. 
Other Subroutines Used 


FUNCTN 


SUBROUTINE RELACO 
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SUBROUTINE EELACO 


General Description 

KEIACO is a FORTRAN subroutine which returns the time required by a 
relay type CO of Westinghouse to close a tripping circuit when the 
st^tings of the relay and the current are specified. 


Calling Sequence 

CALL RELACO(XR,NTCHR,TDS,YR) 

where 

XR is the value of the current expressed in multiples of the pickup 

value. 


NTCHR. is a code for the fvpe of the relay and asstimes integer values 
as follows : 


NTCHR = 1 
2 

3 

4 

5 

6 
7 


type CO- 2 
type CO-5 
type CO- 6 
type CO- 7 
type CO- 8 
type CO- 9 
type CO-11 


TPS is the dial setting and assumes real values in the 


field (0.5 - 11.0) 


YR is the variable which will assume the value of the time to 
close the tripping circuit in the transformed system of coordinates. 
Therefore, the same time in seconds will be: 


time to close the tripping circuit =.01*(10.0**(5,0*YR)) 
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Comments 

For each type of curves, (for example curves CO-5), the curves 
corresponding to the time dial settings .5, 4.0, and 11.0 have been 
stored (with 13 points each) under the symbolic names (XC,YA), (XC,YB), 
and (XC,yc) respectively. 

Specifically: 

(XG(K),YA(D), TDS = .5 

(XC(K),YB(D), IDS = 4.0 K = 1,13 

(XC(K),YC(I)), TDS = 11.0 

I = 1,13 = type CO-2 

I = 14,26 = type CO-5 
I = 27,39 = type CO-6 
I = 40,52 - type CO-7 
I = 53-65 = type CO- 8 
I = 66-78 = type CO-9 
I = 79-91 = type, CO- 11 

Curves for other time dial settings are approxiwted by interpolation. 
Other Subroutine Used 


FUNCTN 


SUBROUTINE FUSECL 
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SUBROUTINE FUSEGL 
General Description 

FUSEGL is a FORTRAM subroutine which returns uhe total clearing 
time and the minimum melting time of any fuse type CL-13 current- limiting, 
14.4 Kv by I-T-E Imperial Corporation, given the current. 

Calling Sequence 

CATL FUSEGL (XPUSE, YFUSEU, YFUSEL, NFUSE) 
where _ 

XFUSE is the current in amperes 

YFUSEU is the variable which will assume the total clearing time of 
the fuse in the transformed system of coordinates. Therefore, the total 
clearing time in seconds will be: 

Total clearing time (sec) = . 01*(10.**(5,*YFUSEU)) 

YFUSEL is the variable which will assume the minimum melting time 
of the fuse in the transformed system of coordinates. Therefore, the 
minimum melting time in seconds will be: 

Minimum melting time (sec) = ,01*(10.**(5,0*YFUSEL)) 

NFUSE is the code for the particular fuse we are interested in. 

The following table explains the selection of this entry: 

Type CL-13 current -limiting fuse units voltage rating 14.4 Kv 
Reference: I-T-E Imperial Corporation 


No. CL- 13-4 

Date 

Rev. 3 


Feb. 20, 1969 



NFUSE 

B 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

L13 

FUSE SIZE 

1 

j lOE 

15E 

20E 

25E 

30E 

40E 

50E 

65E 

80E 

lOOE 

125E 
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Commeixts 

Total clearing time-current characteristic curves are stored with 8 
points for each curve. The minimum melting time is determined by taking 
into account a -15X tolerance in terms of the current of the stored 
curves. That is to find the minimum melting time-current curve, the total 
clearing time-current curve is shifted to the left by a distance corres- 
ponding to 15 percent. 

Other Subroutines Used 


FUNCTN 


SUBROUTINE FUSCLE 
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SUBROUTINE FUSCLE 
General Description 

FUSCLE is a FORTRAN subroutine which returns the total clearing 
time and the minimum melting time of any fuse type CLE-1 and CLE -23 14.4 Kv 
by Westinghouse Electric Corporation, given the current. 

Calling Sequence 

CALL FUSCLE (XFUSE, YFUSEU, YFUSEL, NFUSE) 

where 

XFUSE is the current in amperes 

YFUSEU is the variable which will assume the total clearing time 
of the fuse in the transformed system of coordinates. Therefore, the 
total clearing time in seconds will be 

Total clearing time (sec) = •01'*(10._**(5.*YFUSEU)) 

YFUSEL is the variable which will assume the minimum melting time 
of the fuse in the transformed system of coordinates. Therefore, the 
minimum melting time in seconds will be 

Minimum melting time (sec) =.0l* (10. (5 .0*YFUSEL) ) 

NFUSE is the code for the particular fuse we are interested in. The 
following table explains the selection of this entry 

Type CLE-1 and CLE-2 current limiting power fuses, 14.4 Kv by 
Westinghouse 

References: Curve No. 14 

and Curve No. 15 

Reference No. 598783, 598784, 622124, 622129 




2 

3 

n 

5 

6 

7 

8 

9 

FUSE SIZE 

30E 



40E 



50E 

65E 

80S 

lOOE 

125X 

150E 

175E/200X 
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Comments 

Both minimum, melting time-current and. total clearing time-current 
characteristic curves are stored with 13 points for each curve. In cal- 
culating the minimum melting time, a safety band is taken into account. 
It is achieved by shifting to the right the value of the current by an 
appropriate amount (.05 in the transformed current axis). 

Other Subroutines Used 
PUWCTN 




SUBROUTINE 


FUSESM 
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SUBROUTINE FUSESM 
General Descriptioii 

FUSESM is a FORTRAN subroutine which returns the total clearing 
time and the rainimum melting time of any fuse type SM by S & C Electric 
Company, given the current. 

Calling Sequence 

CALL FUSESM(XFUSE,YFUSEU,YFUSEL,NFUSE) 

where 

XFUSE is the current amperes 

YFUSEU is the variable which will assume the total clearing time of 
the fuse in the transformed system of coordinates. Therefore, the total 
clearing time in seconds will be 

Total clearing time = ,01* (10. ** (5 . 0*YFUSEU) ) 

YFU5EL is the variable which will assume the minimum melting time 
of the fuse in the transformed system of coordinates. Therefore, the 
minimum melting time in seconds will be 

Minimum melting time = . 01* (10.**(5.0*YFUSEL) ) 

NFUSE is the code for the particular fuse we -re interested in. 

The following table explains the selection of this entry 

I, SM Refill units - SLOW SPEED 
Reference: S & C publication: 

Page 1 of 1 

ITGCMUMBER 119-4 

July 19, 1969 


NFUSE assumes values 1 through 19 as follows; 
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NFUSE 

1 

2 

3 


5 

6 

n 

8 

mm 


FUSE SIZE 

15E 

20E 

25E 

30E 

40E 

50E 

65E 

80E 

lOOE 

125E 


NFUSE 

11 

12 

13 

14 

15 

16 i 

17 

1 

18 

19 _ 

FUSE SIZE 

150E 

175E 

200E 

250E 

300E 

400E 

2-250E 

2-300E 

^ 2-400E 


II. SM Refill units - STANDARD SPEED 

Reference: S & C publication: 

I Page 1 of 1 

I II- - I 

TCC NUMBER 153-4 

July 21, 1969 


NFUSE assumes values 20 through 43 as follows: 


NFUSE 



22 

23 

24 

25 

26 

27 

28 

29 

— 

30 

31 

FUSE SIZE 

3E 

5E 

7E 

lOE 

13 E 

15E 


25E 

30E 

40E 

50E 

65E 


32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

80E 

lOOE 

125E 

150E 

17 5E 

200E 

25 OE 

300E 

400E 



2-400E 


Coniment s 

Only the minlmiim melting time - current characteristic curves of 
the fuses are stored with ten points for each curve. Maximum clearing 
time - current characteristic curves are determined by shifting the 
stored curves parallel to the current axis at a distance equal to the 
manufacturer's specified tolerance. 

Other Subroutines Used 


FUNCTN 
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SUBROUTINE FUSEE J 
General Descri-ption 

FUSEEJ is a FORTRAN subroutine which returns the total clearing time 
and the minimuni melting time of any fuse type EJ by General Electric, 
given the current. 

Calling Sequence 

CALL FUSEEJ (XFUSE, YFUSEU, YFUSEL, NFUSE) 

where 

XFUSE is the current in amperes 

YFUSEU is the variable which will assume the total clearing time of 
the fuse in the transformed system of coordinates. Therefore, the total 
clearing time in seconds will be 

Total clearing time (sec) = .01*(10.**(5.*YFUSEU)) 

YFUSEL is the variable which will assume the minimum melting time 
of the fuse in the transformed system of coordinates. Therefore, the 
minimum melting time in seconds will be 

Minimum melting time (sec) = . 01*‘(10.**(5 . 0*YFUSEL) ) 

HFUSE is the code for the particular fuse we are interested in. 

The following table explains the selection of this entry: 

I. CURRENT - LIMITING POWER FUSE 

EJ-1 & EJO-1 14.4 Kv 

References: GES-8104 

GES-8105 

NFUSE assumes values 1 through 9 as follows: 


m 
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II. CURRENT - LIMITING POWER FUSE 
EJ-2 2.4 & 4.8 Kv 

References: GE3-8100A 

GES-8101A 

NFUSE assumes values 20 through 29 as follows; 


NFUSE 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

FUSE SIZE 

2R 



3R 

4R 

6R 

9R 

12R 

18R 

24R 

30R 

36R 


Comment s 

Both miniitium melting time-current and total clearing time-current 
characteristic curves are stored with eleven points for each curve. 

Other Subroutines Used 


FUNCTN 
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SUBROUTINE BKKRDB 
General Description 

BSKRDB is a FORTRA-H subroutine which returns the masimum tripping 
time and the mininatHi tripping time of any breaker type DB by "Westinghouse, 
given the current as a multiple of the pickup current. The subroutine 
does not include special features as ‘'long delay 25-150 seconds,” 

Calling Sequence 

GAEL BRKRDB (NSTD, NLTDB, NSTDE, XLTGR, XSTS, XBR, YUPFER, YLOWER) 

where 

KSTD is the code for the particular breaker we are interested in, 
NSTD assumes values 1, 2, 3 as follows: 


NSTD 

1 

2 

3 

BREAKER 

DB-15 or DB-25 

DB-50 

DB-75.or DB-100 


NLTDB is the long delay setting which should assume integer values 
in the interval (20, 30). Once again, the subroutine features only long 
delay setting adjustable with calibrated marks at 30 and 20 seconds. 

NSTDB is the short delay setting adjustable with calibrated marks 
at 30, 14, and 6 cycles or instantaneous. It assumes integer values 0, 
1, 2, 3 as follows: 


NSTDB 

0 

1 

2 

3 

Short delay setting 

Instantaneous 

6 cycles 

14 cycles 

30 cycles 
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XLTCR is the long delay pickup setting adjustable with calibrated 
marks at 80, 100, 120, 140, and 160 percent of trip unit rating. It 
assumes real values .8, 1., 1.2, 1.4, and 1.6 as follows: 


XLTCR 

.80 

1.00 

1.20 

1.40 

1.60 

Long delay 
pickup setting 

80% 

100% 

120% 

140% 

160% 


XSTS is the instantaneous pickup setting or the short delay pickup 
setting as a multiple of the trip unit rating. This entry assumes a real 
value, 

XER is the current in multiples of the trip unit rating. It is a 
real number. 

YUPPER is the variable which will assume the maximum tripping time 
in the transformed system of coordinates. Therefore, the maximum tripping 
time in seconds will be: 

Maximum tripping time (sec) = .0l*(10.0**(5.0*YUPPER)) 

YLOWER is the variable which will assume the minimum tripping time 
in the transformed system of coordinates. Therefore, the minimum tripping 
time in seconds will be: 

Minimum tripping time (sec) = .01*(10.0**(5.0*YIiOWER)) 

Comments 

Data are stored in the subroutine as follows: 

Total clearing time Minimum tripping time 

(XU (I), YU(D) XL(I), YL(I) 

1=1,9 DB-15 or DB-25 

long delay at 100% pickup and 30 seconds long delay setting 

1=10,15 DB-15 or DB-25 

short delay 6 cycles 
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/ 


1=16,21 DB-15 or DB-25 

short delay 14 cycles 

1=22,27 DB-15 or DB-25 

short delay 30 cycles 

1=28,36 DB-50 

long delay at 100% pickup and 30 seconds long delay setting 

1=37,42 DB-50 

short delay 6 cycles 

1=43,48 DB-50 

short delay 14 cycles 

1=49,54 DB-50 

short delay 30 cycles 

1=55,63 DB-75 or BB-100 

long delay at 100% pickup and 30 seconds long delay setting 

1=64,69 DB-75 or DB-100 

short delay 6 cycles 

1=70,75 DB-75 or DB-100 

short delay 14 cycles 

1=76,81 DB-75 or DB-100 

short delay 30 cycles 


Instantaneous 
Total clearing time 
(XIU(I), YIU(I)) 

1=1,7 DB-15 or DB-25 

1=8,14 DB-50 

1=15,21 DB-75 or DB-100 

Other Subroutines Used 




FUNCTW 


SUBROUTINE BRKEAK 
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SUBfiOUTIWE BRKRAK 
General Description 

BRKSAK is a FORTKAH subroutine which returns the maxiimim and the 
minimum tripping time of any breaker type AK equipped with a series trip 
device EC-IB or EC-I given the current as a multiple of the rated pickup 
current . 

Calling Sequence 

CALL BRKRAK(NSTD, WLTDB, NSTDB, X.,.. ’,R, XBR, YUPPER, YLOWER) 

where 

MSTD is the code for the particular series trip device we are 
interested in, NSTD assumes integer values 1 or 2 as follows: 


NSTD 

1 

2 

Series trip device 

EC-1 

EC- IB 


MLTDB is the code for the long delay setting which can assume 
integer values 1, 2, or 3 as follows: 


NLTDB 

1 

2 

3 

Long time delay band 

Minimum 

Intermediate 

Maximum 


NSTDB is the short time delay setting which assumes integer 
values 1, 2, or 3 as follows: 


NSTDB 

1 

2 

3 

Short time delay band 

Minimum 

Intemediate 

Maximum 


XLTCR is the long time coil rating in per unit of the rated pickup 


current. This entry assumes a real value. 
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XSTS is the short time delay pickup setting expressed in multiples 
of the trip unit rating. This entry assumes a real value 

YUPPER is the variable which will assume the maximum tripping time 
in the transformed system of coordinates. Therefore, the maximum tripping 
time in seconds will be: 

Maximum tripping time (sec) = .01*(10. 0**(5.0*yUPPER) ) 

YLOMER is the variable which will assume the minimum tripping time 
in the transformed system of coordinates. Therefore, the minimum tripping 
time in seconds will be: 

Minimum tripping time (sec) — ,01*(10. 0**(5.0*YLOWER) ) 

Comments 

Data stored in the subroutine are as follows: 

Total clearing time 
XB(I), YB(I) 

Minimum tripping time 
XA(I), YA(I) 


where 


I = 1,13 series trip device EC-1 

@ Minimum long time delay 
Minimum short time delay 

I = 14,26 series trip device EC-1 

@ Intermediate long time delay 
(? Intermediate short time delay 


I = 27,39 series trip device EC-1 

@ Maximum long time delay 
@ Maximum short time delay. 


I = 40,52 series trip device EC-lB 
(? Minimum long time delay 
(3 Minimum short time delay 

\ - 53,65 series trip device EC- IB 

@ Intermediate long time delay 
(3 Intemediate short time delay 

I = 66,78 series trip device EC-lB 
(3 Maximum long time delay 
(3 Maximum short time delay 


Other Subroutines Used 


FUNCTN 


SUB ROUT IHE BRICRLA 
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SUBBDUTINE BRKRLA 
General Description 

BRKRLA is a FORTRAN subroutine which returns the maximum and the 
minimum tripping time of any breaker type LA equipped with a trip device 
Model D or Model A by Allis-Chalmers , given the current as a multiple of 
the pickup current. 

Calling Sequence 

CALL BRKRLA (WSTD, NLTOB , NSTDB, XLTCR, XSTS, YUPPER, YLOWER) 

where 

NSTD is the code for the particular trip device we are interested in. 
HSTD assumes integer values 1 or 2 as follows : 


NSTD 

1 

2 


Model D or DC 

Model 'A or AG 


NLTDB is the long delay setting which should assume integer values 
Ij 2, or 3 as follows: 


NLTDB 

1 

2 

3 

Long time band 

Minimum 


Maximum 


NSTDB is the short delay setting which assumes integer values 1, 2, 
or 3 as folloi'/s: 


NLTDB 

1 

2 

3 

Short time band 

Minimum 

Intermediate 

Maximum 


This entrj' applies only to the Model D, In case of the Model A, this 
entry can assume any integer value without affecting anything- 
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XLTCR does not apply in the present subroutine. Set XLTCR= 1.0. 

XSTS is the instantaneous pickup setting or the short delay pickup 
setting as a multiple of the trip unit rating. This entry assumes a 
real value. 

XBR is the current in multiples of the trip unit rating. It is a 
real number. 

YUPPER is the variable which will assume the maximum tripping time 
in the transformed system of coordinates. Therefore, the maximum tripping 
time in seconds will be: 

Maximum tripping time (sec) = .01*(10.0**(5.0*yupPER)) 

YLOWER is the variable which will assume the minimum tripping time 
in the transformed system of coordinates. Therefore, the minimum tripping 
time in seconds will be: 

Minimum tripping time (sec) = . 01*(10. 0**(5 . 0*YLOWER) ) 

Coinments 

Data stored in the subroutine are as follows: 


Total clearing time 
(XU(K), YU(D), K=l,9 
Minimum tripping time 
(XL(L), YL(M)), L=l,10 


where 


1 = 19 

„ Model D, minimum long time band 

N = 1,10 

N = 11*20 intermediate long time band 


I = 19,27 
N = 21,30 


Model D, maximum Long time band 


I 

N 


28,36 

31,40 


Model A, minimum long time band 


I 

N 


37,45 

41,50 


Model A, intermediate long time band 


I 

N 


46,54 

51,60 


Model A, maximum long time band 


Total clearing time 


XIU(I), YIU(I) 

I = 1,5 Model A, 


YIDL(l) 

YIDU(l) 


Model D, 


YIDL(2) 

YIDU(2) 


Model D, 


YIDL(3) 

YIDU(3) 


Model D, 


instantaneous tripping 
minimum short time band 


intermediate short time band 


maximum short time band 


Other Subroutines Used 


FUNCTN 


SUBROUTINE MCCBGE 
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SUBROUTINE MCCBGE 
General Description 

MCCBGE is a FORTRAN subroutine which returns the maximum tripping 
time and the minimum tripping time of Molded Case Circuit Breakers by 
General Electric (namely ElOO line, F225 line, J600 line and K1200 line), 
given the current as a multiple of the pickup current. 

Calling Sequence 

CALL MCCBGE (NSTD, KLTDB, NSTDB, XLTCR, XSTS, XBR, YUPPER, YLOWER) 

where 

NSTD is the code of the particular breaker we are interested in. 

NSTD assumes values 1 through 14. The selection of this entry is 
explained in note 1. 

NLTDB is the entry which refers to the temperature compensation. 

It assumes the value 1 if the breaker is enclosure compensated; 2 if the 
breaker is ambient compensated. 

NSTDB This entry does not apply to the present subroutine. It is 
kept, however, in order to have the same format for all the breaker 
subroutines. Give the value 1 to this entry. 

XLTCR is the ambient temperature in (°C) . In cases where there is 
difficulty in estimating the ambient temperature, set XLTCR equal to 40.0. 

XSTS is the setting for the instantaneous trip unit and applies in 
the case of the breakers F225, J600, and K1200 only. The possible values 
of this entry appear in note 1. 

In case of breaker ElOO, set XSTS equal to 1. 

XBR is the current in multiples of the pickup current . 
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YUPPER is the varl^le which will assume the maximum tripping 
time in the transformed system of coordinates. Therefore, the maximum 
tripping time in seconds will he: 

Maximimi tripping time (sec) = . 01*(10 ,0**(5. 0*YUPPER) ) 

YLOWER is the variable which will assume the minimum tripping time 
in the transformed system of coordinates. Therefore, the minimum tripping 
time in seconds will be: 

Minimum tripping time (sec) ~ ,01*(10.0**(5.0*YLOWER)) 

Kote 1; 

I. ElOO line by General Electric 

The entry NSTD assimies the following values: 

Types TEB, TED-, and TEED 

NSTD = 1 when the current rating is 15-45 amperes 

= 2 when the current rating is 50-80 amperes 

= 3 when the current rating is 90 5c 100 amperes 

Types TE? 5c THEE (enclosure compensated) 

NSTD = 4 when the current rating is 10-45 amperes 

= 5 when the current rating is 50-100 amperes 

Types TEF 5c THEE (ambient compensated) 

NSTD = 6 when the current rating is 15-40 amperes 

= 7 when the current rating is 50-100 amperes 

II. F225 line, J600 line, K1200 line by General Electric 
The entry NSTD assumes values as follows: 


NSTD 

current rating 
(in amperes) 

possible values of 
the entry XSTS* 

8 

70, 80 

8.5 thru 13.0 


100, no 

6.0 thru 13.0 

9 

90 

6,5 thru 10.0 

125 

5.0 thru 10.0 


150,175,200,225 

4.5 thru 10.0 

J6 

DO 

HSTD 

current rating 
(in amperes) 

possible values of 
the entry XSTS* 

10 

125-400 

3.0 thru 10.0 

11 

450-600 

3.0 thru 10.0 


NSTD 

current rating 
(in amperes) 

possible values of 
the entry XSTS* 

U 

300,350,400 

450,500,600 

3.0 thru 10.0 

13 

700 

3.0 thru 10. 0 


800 

3.0 thru 9.0 

14 

1000 

1 

3.0 thru 7,5 


1200 

3 . 0 thru 6,0 


* The maxinnira value corresponds to the magnetic setting HI and the minimum 
value corresponds to the magnetic setting LO. 


Comments 

Breakers ElOO line: 

Both minimum tripping time-current and maximum tripping time-current 
characteristic curves are stored with 9 and 11 points respectively under 
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the symbolic name (YL, XL) for the minimum tripping time-current curve 
and (YU, XU) for the maximum tripping time-current curve. 

Breakers F225, J6Q0. and KlOO line: 

Long time delay curves and instantaneous curves have been stored 
separately, under the symbolic names [(XAL, YAL) , (XAU, YAU)] and (XIU, YIU) 
respectively . 

Other Subroutines Used 


FUKCm 


STJBR0UTIJ5E BEKROD 
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SUBROUTINE EHKROD 
General Description 

BKKROD is a FORTRAN subroutine which returns the ma: 'mna tripping 
time and the minmum tripping time of any breaker type OD by I-T-E, given 
the current as a multiple of the pickup current. 

Galling Sequence 

CALL BBKROD (NSTD, NLTDB, NSTDB, XLTCR, XSTS,, XBR, YUPFER, YLOWER) 

where 

NSTD is the code for the particular breaker we are interested in. 
KSTD assumes integer values 1, 2 j or 3 as follows; 


NSTD 

1 

2 

3 

BREAKER 

0D~3 & OD-300 

OD-4 & OD-400 

OD-5 & OD-500 
OD-6 6c OD-600 


NLTDB is the long delay setting which assumes integer values 1, 2, 
or 3 as follows; 


NLTDB 

1 

2 

3 

Long-time 
Delay Band 

Minimum 

Intermediate 

Maximum 


In case the manufacturer does not specify the long-time delay band 
(types OD-3 and OD-300), set NLTPB = 1, 

MSTDB is the code for the short time band and assumes integer values 
0, 1, 2 , or 3 as follows: 


NSTDB 

0 

1 

2 

3 

Short-time 

band 

Instantaneous 

Minimum 

Intermediate 

Maximum 
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XLTCR is the long time coil rating in per unit basis. Suggested 
values of this entry ~ .8,, 1.0, 1.2, 1.4, or 1.6. Note that the factory 
setting is XLTCR=1.0. 

XSTS is the instantaneous pickup setting or the short delay pickup 
setting as a multiple of the trip unit rating. This entry assumes a real 
value. 

XBR is the current in multiples of the trip unit rating. It is a 
real number. 

YUFPER is the variable which will assume the maximum tripping time 
in the transformed system of coordinates. Therefore, the maximum tripping 
time in seconds will be: 

Maximum tripping time (sec) = .01*C10.0**(5.0*Y1IPPER)) 

YLQWER is the variable which will assume the minimum tripping time 
in the transformed system of coordinates. Therefore, the minimum tripping 
time in seconds will be: 

Minimum tripping time (sec) - .01*(10.0**(5.0 *YTjOWER)) 

Comments 

Data stored i.; the subroutine are as follows: 

Maximum tripping time 
XLU(I), YLU(I) 

Minimum tripping time 
XLL(I), YLL(I) 

1-1,8 OD-3 and OD-300, long time band 

Same as above. Their existence eliminates a number of 
statements . 

I = 25,32 OD-4 and OD-400, minimum long time delay band 


I = 9,16 1 
I = 17,24J 
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I = 33,40 

OD-4 

and 

I = 41,48 

OD-4 

and 

I = 49,56 

OD-5 

and 

0D“6 

and 

I = 57,65 

OD-5 

and 

OD-6 

and 

I = 65,72 

OD-5 

and 

OD-6 

and 


OD-400, intermediate 
OD-400, maxintum long 

OD-SOO^^, . ^ 

OD- 600 -^^'-’^^™ 

00-500 1^ . ^ 

0D-500\„, . , 

OD-600-'“““ ^“8 


long time delay band 
time delay band 

time delay band 
long time delay band 
time delay band 


Maximum clearing time 
XIU(I), YIU(I) 

I = 1,4 Instantaneous direct acting trip device (curve for XSTS=5.0). 
All types. 


Maximum tripping time 
XSU(I), YSU(I) 


Minimum tripping time 
XSL(I), YSL(I) J 


Short time delay pickup setting: 


XSTS =4,0 


I = 1,6 
I = 7,12 
I = 13,13 


I 

Q 

O 

and 

OD-400 

OD-5 

and 

OD-500 

OD-4 

and 

OD-400 

OD-5 

and 

OD-500 

OD-4 

and 

OD-400 

OD-5 

and 

OD-5 00 


} 

} 

} 


Minimum short time band 


Intermediate short time band 


Maximum short time band 


At 


Other Subroutines Used 


FUNCTN 
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APPENDIX 2 

COORDINATION INDEX COMPUTATION SUBROUTINES 

In this section the subroutines RELFUS, RELBRK, FUSBRK and BKKBRK 
are included. They calculate the coordination indices A2, A1 and RSC, 
as has been described in part 2, for the pairs of devices, relay- fuse, 
relay-breaker, fuse-breaker and breaker-breaker respectively, given the 
type and settings of each device. 

Their numerical valuation of the coordination quality between two 
devices is very useful for the optimization algorithm. 

Detailed description of these subroutines follows: 


SUBROUTINE FUSBRK 


85 


SUBRQUTIjNE fusbrk 
G eneral Description 

FUSBRK is a FORTRAN subroutine v/hich calculates the coordiation 
indices Al, A2, and RSC between a fuse and a breaker when the fuse and 
the settings of the breaker are specified. 

Calling Sequence 

CALL FUSBRK (SHORT, RATIO , ITYPFS , NFUSE , ITYPBR,NSTU ,IiLTDB ,NSTDB ,XLTCR, 
XSTS, PICKUP) where 

SHORT is the short circuit capacity at the point of application of 
the breaker or the interrupting capacity of the breaker whichever is smaller. 

RATIO is the ratio of the current through the breaker and the 
current through the fuse under normal operating conditions. 

The following entries of the subroutine refer to the fuse 

ITYPFS is the code for the type of the fuse. It assumes integer 
values as follows: 


ITYPFS 

i 

2 

3 

4 

type of fuse 

EJ 

CLE 

CL 

SM 

name of the 

corresponding 

subroutine 

FUSEEJ 

FUSCLE 

FUSECL 

L_ 

FUSESM 


NFUSE is the code for the size of the fuse. The selection of this 
entry is explained in the appropriate subroutine. 

The following entries of the subroutine refer to the breaker. 

ITYPBR is the code for the type of the breaker. It assumes integer 


values as follows; 
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ITYPBR 

1 

2 

3 

4 

5 

type of breaker 

AK 

LA 

DB 

Molded 
case by 
G.E. 

OD 

symbolic name of 
the corresponding 
subroutine 

BEKRAK 

BRKRlA 

BRKRDB 

MCCBGE 

BRKROD 


I 

NSTD is the coda Cor the tripping dcTfice . of the breaker (see 
appropriate subroutine). 

NLTDB is the code for the long time delay band (see appropriate sub- 
routine) , 

NSTDB is the code for the short time delay band (see appropriate 
subroutine) . 

XLTCR is a real number referring to the long time coil rating or 
the ambient temperature and should be selected as explained in the appro- 
priate subroutine. 

XSTS is the short time pickup setting (see appropriate subroutine). 
PICKUP is the rated pickup value of the tripping device expressed 
in amperes . 

The subroutine will return in a common block the following numbers; 
Al: a statistical measure of the relative position of the following 

two curves: (i) minimum clearing time of the fuse and (2) maximum 

tripping time of the breaker. 

A2; expected reliability of coordination over the possible values 
of fault current. 

RSC: reliability of coordination at the short circuit capacity of 

the system or at the interrupting capacity of the breaker whichever 


is smaller. 


















Other Subroutines Used 


FUSEE J 
FUSCLE 
FUSECL 
FUSESM 
BRKRAK 
3RKBIA 
BEKRDB 
MCGBGE 
BRKROD 
FUNCTN 

Connnents 

The computatiGn of the coordination indices is performed under the 
following assumptions. 

(1) The field of the possible fault currents is defined 
as all the currents between the pickup current times 
1,Z of the breaker and the short circuit capacity of 
the system at the position of the breaker or the inter- 
rupting capacity of the breeker whichever is smaller . 

(2) The field of the possible fault currents is divided 
into NN equal segments (transformed current axis, i.e. 
log scale) and the expectation integral is calculated 

by summing up the normalized reliability of coordination 
on each segment. 

Suggested value of NN is 20. 

(3) The time to operate of both devices is f-cnsidered 
equally likely distributed between the maximum and 

the minimum clearir time (on the transformed time axis). 


SUBROUTINE BEKBRK 


89 


SUBROUTINE BRKBRK 
General Description 

BRKBRK is a FORTRAN subroutine which calculates the coordination 
indices Al, A2, RSC between two breakers when the settings of the breakers 
are specified. 

Galling Sequence 

CALL BRKBEK(SHORT , RATIO, ITPBRl ,NLTDB1 ,NSTDB1 ,XLTCR1 ,XSTS1 , PCKUPl , 
ITPBR2 ,NSTD2 ,NLTDB2 ,NSTDB2,XLTCR2 ,XSTS2 ,PCKUP2) 

Note: We call breaker #2 the breaker which should trip first and 

breaker #1 the other breaker. Symbolic names referring to 
breaker #2 end with ’’2" and symbolic names referring to 
breaker #1 end with *'l". 

SHORT is the short circuit capacity at the point of application of 
the breaker #2, or the interrupting capacity of either breaker whichever 
is smaller. 

RATIO is the ratio of the current through the breaker #2 and the 
current through the breaker #1 under normal operating conditions. 

The following entries of the subroutine refer to breaker #1. 

ITPBRl is the code for the type of the breaker /rl. It assumes integer 
values as follows: 


ITPBRl 

1 

2 

3 

4 

5 

type of breaker 

AiC 

LA 

DB 

Molded 
case by 

G.E. 

OB 

symbolic name of 
the corresponding 
subroutine 

BRKRAK 

BRKRIA 

BRKRDB 

MCCBGE 

ERKROD 
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NSTDl is the code for the tripping device of the breaker #1 (see 
appropriate subroutine). 

NLTDBl is the code for the long time delay band of breaker #1 (see 
appropriate subroutine. 

MSTDBl is the code for the short time delay band of breaker #1 (see 
appropriate subroutine) . 

XLTGRl is a real number referring to the long time coil rating of 
breaker #1 or the ambient temperature and should be selected as explained 
in the appropriate subroutine. 

XSTSl is the short time pickup setting of breaker #1 (see appropriate 
subroutine) . 

FCKUPl is the rated pickup value of the tripping device, breaker #1, 
expressed in amperes. 

The following entries of the subroutine refer to breaker #2. 

Explanation same as for breaker #1. 

ITPBR2 

NSTD2 

NLTDB2 

NSTDB2 

XLTCR2 

XSTS2 

PGKUP2 

The subroutine will return in a common block the following numbers: 
Al ; a statistical measure of the relative position of the following 
two curves: (1) minimum tripping time of breaker #1, and (2) maxi- 

mum tripping time of breaker #2. 

A2: expected reliability of coordination over the possible values 


of fault current. 
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RSG : reliability of coordination at the short circuit cap<acity 

of the system or at the interrupting capacity of either breaker 
whichever is smaller. 

Other Subroutines Used 

BRKEAK 

BRKRLA 

BRKRDB 

MCCBGE 

BSKROD 

FUNCTN 

Comments 

The computation of the coordination indices is performed under the 
following assumptions; 

(1) The field of the possible fault currents is defined as all the 
currents between the maximum pickup current times 1.2 of the 
two devices (i.e. the maximum of the two numbers (1.2)x(r-'=*ted 
pickup of breaker l)x(long time coil rating in per unit of 
breaker 1) and (1.2)x(rated pickup of breaker 2)x(long time coil 
rating in per unit of breaker 2)) and the short circuit capacity 
of the system. 

C2) The field of the possible fault currents is divided into NN 
equal segments (transformed current axis, i.e. log scale) and 
the expectation integral is calculated by summing up the 
normalized reliability of coordination on each segment. 

Suggested value of NK is 20. 

(3) The time to trip of both devices is considered equally likely 
distributed between the maximum and the minimum tripping time. 

(on the transformed time axis). 


SUBROUTINE EELFUS 
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SUBROUTINE EELFUS 
General Descrjptl'i n 

RELFUS is a FORTRAN subroutine which calculates the coordination 
indices Al, A2 and RSC between a relay and a fuse when the fuse and the 
relay setting are specified. 

Cal ling Sequence 

CALL EELFUS (SHORT, RATIO, ITyPRL,NTCHR,TDS,PCKPEL,ITYFFS,NFUSE) 

where 

SHORT is the short circuit capacity at the point of application 
of the fuse. 

RATIO is the ratio of the current through the fuse and the current 
through the relay under normal operating conditions. 

The following entries of the subroutine refer to the relay ITYPRL 
is the code for the type of the relay. It assumes integer values as 
follows : 


ITYPRL 

1 

2 

type of relay 

lAC 

GO 

name of the 

corresponding 

subroutine 

RELIAC 

REIACO 


NTGHR is the code of the family of the relay characteristics 
(see appropriate subroutine). 

i^CKPRL is the pickup current of the relay. 

The following entries of the subroutine refer to the fuse 

ITYPFS is the code for the type of the fuse. It assumes integer 


values as follows: 
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ITYPFS 

1 

2 

3 

4 

type of fuse 

EJ 

CLE 

CL 

SM 

name of the 
c orr e s p onding 
subroutine 

FUSEEJ 

FUSCLE 

FUSECL 

FUSESM 


NFUSE is the code for the size of the fuse. The selection of 
this entry is explained in the appropriate subroutine. 

The subroutine will return in a connnon block the following numbers; 
Al : a statistical measure of the relative position of the 

following two curves: (1) mean operating time of the relay 

and (2) maximum clearing time of the fuse, 

A2; expected reliability of coordination over the possible 
values of fault current. 

RSC : reliability of coordination at the short circuit capacity 

of the system. . 


Other Subroutines Used 


RELUC 

RELACO 

FUSEEJ 

FUSGLE 

FUSECL 

FUSSSM 

FUNCTN 


Comments 

The computation of the coordination indices is performed under the 
following assumptions: 

(1) The field o^ the possible fault currents is defined as all the 
currents between the pickup current times 1.2 of the relay and 
the short circuit capacity of the system at the position of the fuse. 
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(2) The field of the possible fault rurreixts is divided into NN 
equal segments (transformed current axis, i.e. log scale) 
and the expectation integral is calculated by summing up the 
normalized reliability of coordination on each segment. 

Suggested value of NN is 20, 

(3) The relay time to trip is considered equally likely distributed 

on the interval t' i .03 where t’ is the nominal value on the 
N N 

transformed time axis. The fuse time to melt is considered 
equally likely distributed between the minimum melting time 
and the maximum clearing time. 




O' 


SUBROUTINE RELBEK 
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SUBROUTINE BELBHK 
General Description 

RELBEK is a FORTRAN subroutine which calculates the coordination 
indices Al, A2 and R'3C between a relay and a breaker when the relay and 
breaker settings tre specified. 

Calling Sequence 

CALL RELBRK (SHORT, RATIO, ITYPRL,NTCHR,TDS,PCKPRL,ITYPBR,NSTD,NLTDB, 
NSTDB ,XLTGR,XSTS ,PCKPBR) 

where 

SHORT is the short circuit capacity at the point of application 
of the breaker or the interrupting capacity of the breaker whichever is 
smaller. 

RATIO is the ratio of the current through the breaker and the 
current through the relay. Under normal operating conditions. 

The following entries of the subroutine refer to the relay, 

ITYPRL is the code for the type of the relay. It assumes integer 
values as follows : 


ITYFRL 

1 

2 

type of relay 

lAC 

CO 

name of the 

corresponding 

subroutine 

RELIAC 

RELAGO 


NTCHR is the code of the family of the relay characteristics (see 
appropriate subroutine) . 

PCKPRL is the pickup current of the relay. 


i Us_ 



The following entries of the subroutine refer to the breaker. 

ITYPBR is the code for the type of the breaker. It assumes integer 
values as follows: 


ITYPBR 

1 

2 

3 

4 

5 

' ' 

type of 
breaker 

AK 

LA 

DB 

folded, 
case by 

G.E. 

OD 

symbolic name of 
the corresponding 
subroutine 

BRKRAK 

BRKRLA 

BBKKDB 

MCCBGE 

BRKROD 


HSTD is the code for the tripping device of the breaker (see 
appropriate subroutine) . 

NLTDB is the code for the long time delay band (see appropriate 
subroutine) . 

NSTDB is the code for the short time delay band (see appropriate 
subroutine) . 

XLTCR is a real number referring to the long time coil rating or 
the ambient temperature and should be selected as explained in the appro- 
priate subroutine, 

XSTS is the short time pickup setting (see appropriate subroutine). 
PCKPBR is the rated pickup current of the breaker. 

The subroutine will return In a common block the following numbers: 
Al : a statistical measure of the relative position of the 

following two curves: (1) mean operating time of the relay and 

(2) maximum tripping time of the breaker. 

A2 : expected reliability of coordination over the possible 

values of fault current. 

RSG : reliability of coordination at the short circuit capacity 

of the system. 
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Other Subroutines Used 

RELIAC 

REIAGO 

BKKRAK 

BRKRIA 

BEKRDB 

MGCB6E 

BRKROD 

FUNCTN 

Comments 

The computation of the coordination, indices is perfomed under the 
following assumj»tions : 

(1) The field of the possible fault currents is defined as all 
the currents between the maximum pickup current times 1.2 of 
the two devices (i.e. the maximum of the two numbers (l,2)x 
(relay pickup current )x (ratio) and (1.2)x (rated pickup of 
the breaker)x(long time coil rating in per unit of the 
breaker)) and the short circuit capacity of the system. 

(2) The field of the possible fault currents Is divided into M 
equal segments (transformed current axis, i.e, log scale) 
and the expectation integral is calculated by summing up 
the normalized reliability of coordination on each segment. 

Suggested value of HN is 20. 

(3) The relay time to trip is considered equally likely distributed 
on the interval t^ i .03 where is the nominal value on 

the transformed time axis. The breaker time to trip is con- 
sidered equally likely distributed between the maximum and 
the minimum tripping time (on the transformed time axis). 
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APPENDIX 3 

OVERHEATING AND MOTOR STARTING CURRENTS 
PROTECTION SUBROUTINES 

Two subroutines are presented here. The first one, CHECKH, 
assesses protection of the system against overheating. It acttially 
checks if the current limiting fuse employed in the system can clear 
before a fault current can cause any heat damage to the system. The 
second one, CHECKS, assesses protection of the system against unnecessary 
tripping of the breakers due to motor starting inrush currents. It 
actually checks if the breaker closest to the motor will not trip under 


normal motor starting conditions. 

Detailed description follows: 


SUBROUTINE CHECKH 


SUBROUTINE CHECKH 


General Description 

CHECKH is a FORTRAN subroutine which checks if the fuse clears in 
less time than that which might damage the system* 

Galling Sequence 

GALL CHECKH (SHORT, ITYPFS,NFUSE,NGODE,NPASS) 

where 

SHORT is the short circuit capacity of the system at the point of 
application of the fuse. 

ITYPFS is the code for the type of the fuse. It assumes integer 
values as follows: 


ITYPFS 

1 

2 

3 

4 

type of fuse 

EJ 

CLE 

CL 

SM 

name of the 

corresponding 

subroutine 

FDSEEJ 

FUSCLE 

FUSEGL 

FUSESM 


NFUSE is the code for the size of the fuse (see appropriate sub- 
routine) . 

NGODE offers an option; 

if NCODE " 0, no check will be made and the fuse will "pass.” 

Any other value of NCODE will yield in checking if the fuse 
meets the protection requirements. 

NPASS is the code which will be returned by the subroutine. 

If NPASlJ equals 2, the fuse does not meet the protection 
requirements. 

If NPASS equals 1, the fuse passes the protection requirements. 
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Comments 

The designer should provide the following information: The heat 

curve defined with no more than 20 points, under the symbolic name (XA,YA). 

XA; current in amperes 
YA: time in seconds 

Other Subroutines Used 

FUNCTN 

FUSSEJ 

FUSCLE 

FUSEGL 

FUSESM 

Note: In almost all cases of the KSG power system the prohibitive heat 

dissipation in case of a fault will occur in the nearest transformer. 
Therefore the designer should determine the heat curve of the transformer 


in most of the cases. 


SUBROUTINE CHECKS 



SUBROUTINE CHECKS 


General Description 

CHECKS is a FORTRAN subroutine which checks i£ the jninimtan tripping 
curve of a feeder breaker is above the starting current curve (time vs 
current) of a motor. 

Calling Sequence 

CALL CHECKS (SHORT, ITYPBR,NSTD,NLTDB,NSTOB,XLTCR,XSTS, PICKUP ,NCODE, 
NPASS) 

where 

SHORT is the short circuit capacity of the system or the interrupting 
capacity of the breaker whichever is smaller. 

The following entries of the subroutine refer to the breaker. 

ITYPBR is the code for the type of the breaker. It assumes integer 
values as follows: 


ITYPBR 

1 

2 

3 

4 

5 

type of 
breaker 

AK 

LA 

DB 

Molded 
case by 

G.E. 

OD 

sjrmbolic name of 
the corresponding 
subroutine 

BSKRAK 

BRKBIA 

BRKRDB 

MCCBGE 

BRKROD 


NSTD is the code for the tripping device of the breaker (see 
appropriate subroutine). 

NLTDB is the code for the long time delay band (see appropriate 
subroutine . 

HSTDB is the code for the short time delay band (see appropriate 


subroutine) . 
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XLTGR is a real number referring to the long time coil rating 
or the ambient temperature and diould be selected as explained in the 
appropriate subroutine. 

XSTS is the short time pickup setting (see appropriate subroutine) , 
PICKUP is the rated pickup value of the tripping device expressed 
in amperes. 

NCQDE offers an option: 

If NCOUE = 0, no check will be made and the breaker will "pass," 

Any other value of this entry will yield a normal run of the 
subroutine. 

NPASS is code which will be returned by the subroutine 

If NPASS = 1 (2), the minimum tripping curve of the breaker 
lies (does not lay) above the starting current durve of the 
motor . 

Comments 

The designer should provide the starting current curve of the motor 

under the symbolic name (XB,YB) and with no more than 20 points. 

Note; XB; current in amperes 
YB: time in seconds 

Other Subroutines Used 

BEKEAK 

BEKELA 

BEKEDB 

MCGBGE 

BEKROD 

FUNCTN 
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APPENDIX 4 
PLOTTING SUBROUTINES 

Two subroutines are included here. The first one can overlay the 
characteristics of at most seven devices, namely, two relays, two fuses 
and three breakers. It has been designed to be used independently. It 
provides an excellent illustration of the coordination status. 

The second subroutine is similar to the one above. It can overlay 
the characteristics of at most four devices, namely, one fuse and three 
breakers. Its format is readily applicable to the most common case of 
substation at the Kennedy Space Center power system. It can be called 
by the optimization subroutine to plot the optimal curves of the employed 
devices (see Appendix 5), 

Detailed discription follows: 


4 


4 


\ 
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k. Plotting Routine 

This routine can plot the characteristics o£ 7 protective devices, 
namely, two relays, two fuses and three breakers. 

It requires 8 data cards. 

The format of the data cards follows; 

1st data card 

1-10 current transformation ratio between 

relay A & B 

11-20 current transformation ratio between 

relay B & fuse A 

21-30 current transformation ratio between 

fuse A & fuse B 

31-40 current transformation ratio between 

fuse B fit breaker A 

41-50 current transformation ratio between 

breaker A fit breaker B 

51-60 current transformation ratio between 

breaker B fit breaker G 

61 code for the location of the substation 

= 1 if in Launch Complex 39 
“ 2 if in Industrial Area 

2nd data card 

Refers to relay A, If relay A does not exist leave the card blank. 
If relay A exists fill in as follows; 

i-5 symbolic name of the bus where the 

relay is applied 

6-15 the relay* s pickup current setting 


16-25 the interrupting capacity of the 

associated breaker 

26-30 the code for the relay tjrpe 

1 for Relay IA.C 

2 for Relay CO 

31-35 the code for the particular 

characteristic of the relay (see 
appropriate subroutine) 

36-40 titfie dial setting of the relay 

3rd data card 

Refers to relay B. Explanation and format are same as for relay A 
4th data card 

Refers to fuse A. If fuse A does not exist leave the card blank. 
If fuse A does exist fill in as follows: 

1-5 symbolic name of the bus where the 

fuse is applied 

6-10 code for the ti^e of the fuse 



= 1 

type 

EJ 


2 

type 

CLE 


3 

type 

CL 


4 

type 

SM 

11-15 

code for 

the size 

of the fuse 

16-25 

current 

rating of 

the fuse 


5th data card 


Refers to fuse B. Format is the same as par fuse A. 


6th data card 


Refers to breaker A. If breaker A does not exist leave the card 
blank. If breaker A doss exist fill in as follows: 


1-5 


6-15 

16-25 

29-30 


31-35 

36-40 

41-45 

46-55 

56-65 


symbolic name of the bus where the 
breaker is applied 

rated pickup setting in amperes 

interrupting capacity of the breaker 
in amperes 

code for the type of the breaker 


^ 1 

type 

AK 

2 

type 

LA 

3 

type 

DB 

4 

type 

Molded case 
by G.E. 

5 

type 

OD 


code for the series trip device 
code for the long time delay band 
code for the short time delay band 
long time coil rating in per unit 
short time pickup current setting 


7th data card 

Refers to breaker B, Format as per breaker A, 


8th data card 


Refers to breaker C. Format as per breaker A 
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The following scheme illustrates the general case covered by 
the plotting routine. 


relay A 


M relay B 


n fuse A 


rn fuse B 


breaker A 


[~~1 breaker B 
i~l breaker G 


Of course, in a specific case some of the referred devices might 
be missing. Insert blank data cards for those devices and the program 
will neglect them. 

Calling Procedure 

To call the plotting routine it is sufficient to provide the eight 
data cards. The device characteristic subroutines should be included 
in the program 


4 




Comments 


For each device, the set of currents applicable to the case is 
defined as all the currents between the pickup value (or rated current 
for fuses) and the short circuit capacity of the system at the position 
of application of the device- This interval is then divided in forty 
equal subintervals (always on the transformed axis, i.e. x = log I+Gp. 
For the middle value of each sub interval the operating time is calculated 
by calling the appropriate device subroutine. Therefore, the curve (or 
the curves - minimum and maximum operating times) of each device is 
determined with forty points each. 

The datagraphics 4460 system is called to put these points on a 
Common coordinate system and then draw a line through the points belonging 
to the same curve. 

Note that the routine is written in FORTRM 


Other Subroutines Used 

RELIAC 

RELACO 

FUSEEJ 

FUSEGL 

FUSCLE 

FUSESM 

BRKEAK 

BRKRIA 

BRKROD 

MCCBGE 

5RKRDB 

FUNGTN 


B. Subroutine BLOT 

Subroutine PLOT is a FOETEAN subroutine which can overlay the time 
vs current characteristics of at most four devices, namely, one fuse, and 


three breakers. 


Calling Sequence 


CALL PLOT (ETFA,RTAB,KTBC,ITPFF,NFUSEO,SHFUSE,PCKPBA,BAINTC, 

ITPBA,NSTDAO ,NLTBAO,NSTBAO ,XLTGAO ,XSTSAO, SHBA, PCKPBB ,BBINTC , 
ITPBB,NSTDBO,HLTBBO,NSTBBO,XLTCBO,XSTSBO, SHBB ,PCKPBG,BCmTG , 
ITPBG jNSTDCO jNLTBCO ,NSTBCO ,XLTCCO , XSTSCO,SHBC) . 

where 

ETFA current transformation ratio between fuse and breaker A 
RTAB current transformation ratio between the breaker A and B 
RISC current transformation ratio between the breakers B and G 
ITPFF code for the type of fuse (see subroutine EELFUS) 

N FUSED code for the size of fuse (see appropriate subroutine) 

SHFUSE the short circuit capacity of the system at the point of 
application of the fuse. 

Following entries refer to breaker A. 

PCKPBA the rated pickup current of breaker A 
BAllCTC the breaker A interrupting capacity 

ITPBA code for the type of breaker A (see subroutine BEKBBK). 

HSTDAO code for the series trip device of breaker A 
NLTBAO code for the long time delay band of breaker A 
USTBAO code for the short time delay hand of breaker A 
XLTCAO the long time coil rating of breaker A (in per unit) 

XSTSAO the short time pickup value of breaker A in multiples of the 
rated pickup circuit 

SHBA the short circuit capacity of the system at the position of 


application of breaker A. 
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Following entries refer to breaker B. Explanation in similar as per 

breaker A* 

FGKPBB 

BBINTC 

ITPBB 

WSTDBO 

5JLTBB0 

NSTBBO 

XLTCBO 

XSTSBO 

SHBB 


Following entries refer to breaker C. Explanation is similar as per 

breaker A 

PCKPBC 

BCIKTC 

ITPBC 

NSTDCO 

NLTBGO 

NSTBCO 

XLTCCO 

XSTSCO 

SHBC 


Coimtients 


This subroutine is part of the plotting program. Explicity, it 
retains the part of the program which refers to fuse B^ breaker A, 
breaker B and breaker C, 


APPENDIX 5 


COORDINATION OPTIMIZATION SUBROUTINES 


General 

The purpose of these subroutines is to determine optimal settings 
of existing systems by utilizing the concepts developed in part 2. 

The accepted procedure is very simple. For an acceptable com- 
bination of devices and setting the coordination indices are calculated. 
Then a weighted average is computed. This average is the basis of com- 
parison. The configuration with the highest weighted average index is 
defined as optimal. 

The average coordination index is defined by the form 
A.C.I. = XA*A2+XB*A1+XC*RSC 
with XA+XB+XC =1.0 

XA, XB and XC can be determined by the user according to his 
philosophy of coordination and according to the nature of the system 
under study. 

For example, if coordination at the short circuit capacity only 

is desired, it is sufficient to set 

XA = 0.0 
XB = 0.0 
XC = 1.0 

The above extreme case reveals that by assigning appropriate 
values to these coefficients it is possible to satisfy different tendencies 
in the coordination philosophy. 
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Main Program FBBBC0 
General Description 

Program FBB0C is written in FORTRAN. It determines the optimma 
coordinated settings of the Involved protective devices by searching 
and comparing all possible settings. The comparison is executed by 
means of a predescribed criterion which can be modified by any user, if 
desirable. In particular the program has been written in connection 
with the following very common case of substations 



Required Data 

Before the data cards the user should provide three cards on 
which he can punch anything he wishes. This text will be the title of 
the output. 

The program requires 7 data cards to be filled in as follows; 

1st data card 

1-6 symbolic name of the bus where the 

fuse is located 

7-12 symbolic name of the bus where the 

breaker A is located 

13-18 symbolic name of the bus where the 

breaker B is located 

19-24 symbolic name of the bus where the 

breaker C is located 


116 


25 

26 

27 

28 

2nd data card 

1-10 

11-20 

21-30 

3rd data card 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


code for option #1 (see options) 

code for option #2 (see options) 

code for the area where the substation 
is located 

= 1 if in Launch Complex 39 
=2 if in Industrial Area 

code for option #3 (see options) 


interrupting capacity 
interrupting capacity 
interrupting capacity 


of breaker A 
of breaker B 
of breaker C 


current transformation ratio between 
fuse and the breaker A 

current transformation ratio between 
the breakers A and B 

current transformation ratio between 
the breakers B and C 

the rated pickup current in amperes 
of breaker A 

the rated pickup current in amperes 
of breaker B 

the rated pickup current in amperes 
of breaker C 
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4th data card 


1-5 

6-10 

11-15 

16-20 


code for the smaller size candidate fuse 

code for the larger size candidate fuse 

step (insert 1 to include all the 
intermediate sizes) 

code for the type of fuse 


1 

for 

fuses 

type 

EJ 

2 

for 

fuses 

type 

CLE 

3 

for 

fuses 

type 

CL 

4 

for 

fuses 

type 

SM 


note: to determine correspondence between fuse size and 

code consult the subroutine which refers to the particular 
type of the fuse. 


5th, 6th, and 7th data cards 

The above cards nominate all the candidate breake?* settings to 
be considered by the program. The format is identical for the three 
cards as follows: 


1-5 

smaller 

6-10 

larger 

11-15 

step 


code for the series trip device 


16-20 

21-25 

26-30 



code for the long time delay band 


31-35 

36-40 

41-45 



code for the short time delay band 
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46-50 

smaller i 

1 


51-55 

larger 

1 the long time coil rating expressed 

in 

56-60 

1 

stap 

J percent (integer) 


61-65 

smaller 

i 

) 


66-70 

larger 

* the short time pickup setting (integer) 

71-75 

step 



76-80 


code for the type of the breaker as 

follows : 


1 for breakers type AK 

2 for breakers t 3 rpe LA 

3 for breakers type DB 

4 for molded case breakers by G.E. 

5 for breaker type OD 

note: the user should study the selection of data for the 

worked example (attached) in order to resolve ambiguities. 

Subroutines Used 


FUSBEK 

BKKBRK 

CHECKS 

CHECKS 

FUSEEJ 

FUSCLE 

FUSEGL 

FUSESM 

BEKSAK 

BEKEIA 

BRKEDB 

MCCBGE 

BEKROD 

FUNCTN 


Options 

option #1 The user can secure protection of the system against overheating 
by inserting any number except zero (0) in column 25 of the first 
data card and providing the necessary data in subroutine CHECKH 
(see subroutine for further information) 
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option #2 If the system includes motors of size commensurate to the 
rating of the system, the program provides the option to secure 
the system against unnecessary tripping due to the motor starting 
currents. In order to do so fill in column 26 of the first data 
card any number except zero (0) and provide the necessary data in 
subroutine CHECKS (see subroutine for further information). 

option #3 This option refers to the type of the substation: 

(1) Double fed substation with primary fuse, main secondary 
breaker, tie line breaker and feeder breakers is recognized with 
the code 0 punched in column 28 of the first data card. 

(2) Single fed substation with primary fuse, main secondary 
breaker and feeder breakers is recognized with the code 1 
punched in column 28 of the first data card. 

In the last case, where the tie line breaker is missing, the sixth 
data card which refers to the missing tie line breaker should be a 
duplicate of the fifth data card. 

Flow Chart 


A flow chart of the program is presented in the next pages. 
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(1) double fed substation 
,(2) single fed substation 




1 




\V 


calculate coordination indices 
between breakers A and B 
(A2, A1 and RSG) 


i BELOC2 = CRTSB 

Set parameters 
of breaker B 
, equal to those 
of breaker A 


b:<-{ 

" i 


i next 
settings 


1 NO 


._1 


.X. 


calculate an overall coordina- i 
tion index between breakers 
A and B 

i’^EL0C2 = XA*X2+XB*X1+XC*RSC 
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Read Data 





„ 

calculate short 
circuit capacities 



1 compare short circuit 
I capabilities with inter 
I rupting capacities of the j 
i br eak ers . | 



— V . 

CRTRA. = .60 
GRTRB = CRTRA 
HELP = CRTRA**3 




y 

assume settings 

.. . * 4 ^,. 


J 


"r 


■>! 


-fi'. 



{ 


next fuse 


calculate coordination 
indices between fuse 
and breaker A 
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check if 

^'^'^breaker C protects'''-.^ 
■^he system against unneces-' 
■'••^sary tripping due to 
' motor starting 
current 


NO 

next ! 


settings | 


YES 



r 


X 


' calculate coordination indices 
j between breakers B and C 
I (Al, A2 and RSC) 


t 


! 


, . 

calculate an overall coordina- 
tion index (between breakers 
B and C) 


REL0C3 XA*A2+XB*A1+XC*RSC 


X 

calculate overall coordination 

index of the system 

RELOC = REL0C1*REL0C2*REL0C3 



I nominate as "updated optimum | 
f settings" the current settings { 
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eeloc 


— 

CRTEA= 


_Y... 


*3 j 

CRTRB = ^LP \ 


L.s^lligs. 


do more 

^settings exist?/ 


,/^is current \ NO 

CTEA greater > 

V than .60? / 


CRTM-CRTRA- . 1 1 
-GRTRB=GRTEA 
RELP=(CRTRA)^ 


YES 



/ print out title, 
optimum settings, short \ 
! circuit capacities, and | 
\ operating times at the / 
\ short circuit capacit^ 
in seconds. 

I 


i 

. ± 

Gall the plotting subroutine 
for the optimal settings of 


1 

i 


the devices 



display of the 
characteristics 
on the Data 
graphics 4460 
..system 


. „V.,. 

STOP. 
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The flow chart shows the mechanics of the program and makes obvious 
that it can be modified to meet other requirements. For example, 
suppose that a relay controlled breaker is in place of the breaker A. 

In this case it would be sufficient to make only two changes : 

(1) Instead of calculating the coordination indices between the 
fuse and the breaker A (subroutine FUSBRK), the coordination 
indices between the fuse and the relay should be computed 

by calling the subroutine KELFTJS. 

(2) Instead of calculating the coordination indices between the 
breakers A and B (subroutine BRKBEK) , the coordination 
indices between the relay and the breaker B should be com- 
puted by calling the subroutine RELBEK. 

With this we encourage users to modify and/or supplement the routines 
according to their needs. 


